Corrections to the Coulomb pseudopotential in the theory of superconductivity 
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We calculate corrections to the conventional result for the Coulomb pseudopotential in the theory 
of superconductivity. We find that retardation effects are still operative for higher order correc- 
tions, but less efficient due to a reduction of the effective bandwidth, leading to larger values of 
the pseudopotential and reduced superconductivity. Our calculation are based on a perturbation 
expansion for the Hubbard-Holstein model, and the conclusions are corroborated by comparison 
with non-perturbative dynamical mean-field results. 
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The theory of superconductivity based on the electron- 
phonon mechanism has been very successful in describing 
the properties of many materials with impressive accu- 
racy [l|, 01 ■ The electron phonon coupling is treated in 
the Migdal Eliashberg (ME) theory [3[, which relies on 
Migdal's theorem Ql- This employs the fact that typi- 
cal electron (Ed) and phonon (wph) energy scales differ 
largely. Then perturbation theory greatly simplifies as 
vertex corrections are small. This is true even for large 
values of the electron-phonon coupling parameter A > 1 
as long as Awph/i^oi remains small [4 6]. 

A crucial issue is the effect of the Coulomb repul- 
sion, since this is typically much larger than the phonon- 
induced attraction. The electronic repulsion in the pair- 
ing channel can be projected to the phonon scale. It is 
then str ong ly reduced due to retardation effects and one 
finds 04io[, 



(1) 



often termed the Morel- Anderson pseudopotential. Here, 
fJ-c ~ PqU, where [/ is a typical screened Coulomb in- 
teraction and pq is the density of states (DOS) at the 
Fermi energy. Since typically E^i 3> i^ph, one finds that 
/i* ^ fie and often also /i* < A. Eq. ([1]) leads to esti- 
mates of the order /i* ^ 0.1 — 0.14. This agrees rather 
well to the fitting parameter fi* obtained from tunneling 
spectroscopy for many conventional superconductors [l| . 

Although ME theory has been very successful, the 
treatment of the Coulomb repulsion is by no means rig- 
orous [lH. For the electron-phonon interaction, Migdal's 
theorem justifies the neglect of vertex corrections. For 
the Coulomb interaction there is no similar justification, 
and Eq. ([T]) is based on an uncontrolled approximation. 

There is experimental evidence that Eq. ([T]) may be 
insufficient for some materials. For, e.g., V or NbsGe [ij, 
the values for /i* appear to be of the order 0.2 — 0.3 sub- 
stantially larger than the traditional quotes, even though 
the ratio Ed/ujph is not much different. Density func- 
tional theory (DFT) finds good agreement with the tun- 
neling results, but to explain the experimental values for 



Tc, quite large values of /i* have to be used [l^- For 
elemental Li at ambient pressure [l3l [13 the coupling 
constant was estimated to be A ~ 0.4 |l3l. With 
fi* ^ 0.1 this implies Tc ~ IK, while experiments recently 
found superconductivity at Tc ^ 0.4mK at ambient pres- 
sure [l^. For alkali-doped fuUerides, the Morel- Anderson 
theory leads to a larg e reduction of p,* due to coup ling s to 
higher subbands 17 1 , although this is unphysical 17 , l3| 
and raises serious questions about the neglect of vertex 
corrections for molecular solids. Actually, it was found 
that superconductivity in fullerides is due to a compli- 
cated interplay between the Coulomb interaction and the 
interaction with Jahn- Teller phonons (l9l. [20|. 

The Morel-Anderson theory corresponds to treating 
the irreducible vertex to first order in U. For the model 
considered below, we extend this to second order in U 
within a projection approach and provide numerical cal- 
culations of fi* without further approximation, as well as 
approximate analytical calculations. We find that retar- 
dation effects lead to the reduction fic ~^ also in the 
higher order calculation, but not as efficiently as for the 
first order one. Then we add the electron-phonon interac- 
tion and calculate the superconducting gap numerically 
and approximately analytically. The results are well un- 
derstood by the derived results for ji* . To see the effect of 
higher order corrections, our calculations are carried out 
such that they can be compared with non-perturbative 
dynamical mean- field theory (DMFT). DMFT often pro- 
vides a good approximation for three dimensional sys- 
tems and fully treats retardation effects, which are cru- 
cial for the Coulomb pseudopotential. Finally, we show 
that results based on the perturbation theory agree well 
with DMFT calculations up to intermediate couplings. 

We choose a minimal model, which possesses all nec- 
essary ingredients, the Hubbard-Holstein model, 

J 



H = -V(t 



+ h.c.) + Uj2 



(2) 



-UJQ 



^5l5,+g^(6, + 5j)(^n,,,-l). 
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c\ creates an electron at lattice site i with spin a, and 
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bl a phonon with oscillator frequency Wq, '^i.o- = cl^Ci^^- 
The electrons interact locally with strength U, and their 
density is coupled to an optical phonon with coupling 
constant g. In the limit of large dimensions this model 
can be solved exactly by the DMFT. 

First we deduce /z* from an analysis of the pairing in- 
stability, in the limit iw„ 0, q 0. We define the 
symmetric matrix [21|, 



0)v^ 



(3) 

where /3 is the inverse temperature, F^pp^ {ii^m , i<jJn2 ] 
is the irreducible vertex in the particle-particle chan- 
nel and the pair propagator reads x'^(iuJn) ~ [G{iujn) — 
G(-iix)„)]/[C(-iw„) - Cii'^n)], where C(«w„) = iw„ + fJ-- 
S(icj„) and G{iujni) is the local lattice Green's function. 
A becomes singular at Tc- We introduce the "low-energy 
part" 



A. 



low 
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■^nn' \_-A \n'm' -^r, 
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n,m such that |a;„|, |a;,„| < Wph, and A is the the block 
for |a;„|, \ujm\ > Wph- If A'"™ is singular, A is also singu- 
lar. The "folding in" of larger frequencies describes how 
retardation effects reduce effects of the Coulomb repul- 
sion on low frequency properties. We first consider the 
lowest order term of F^pp) in U, F^pp'^^ = -U. We focus 
on the dependence on the half-band width D and assume 
a constant DOS, po — 1/{2D). It is a rather good ap- 
proximation to write x'^(iLL>n) = PaT^/\^n\, if l^^nl < D 
and otherwise. With /ic ^ poU, A takes the form. 



nni — ^nm 



l3\/\uJnl. 



(5) 



which is separable and can be inverted exactly. Replacing 
summations by integrals, we find then 



4I0W _ r 



Pc 
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Comparison of Eq. ([5|) and © leads to the Coulomb pseu- 
dopotential, pc ^ Pc, as given in in Eq. ([T|). 

We next consider the second order term of F'pp^ in [/, 
which comes from a crossed diagram. 



where the particle-hole bubble is given by 



(7) 
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We can write n(?'cj„) = —f{x)apQ, x — icUnPo, where / is 
independent of D and approximated as 



1 



1 -f- b\x\ + cx^ ' 



(9) 
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FIG. 1: (Color online) /i* as a function of p^ for D/u)p\i = 100 
and Pu)pYi = 240. The figure shows the calculated results using 
both the first order and first plus second order result for r'''^' 
as well as the approximation in Eq. (|10|) . 



where a = 1.38, b = 2 and c = 5 are suitable values for 
the constant DOS. 

Because of the form of n(icj„ + ia;„i), A in Eq. (|4]) can- 
not be inverted analytically. Instead we use the inverse 
of A based on F'^pp^'^, which is now only correct to first 
order in U. However, since the off-diagonal terms of A 
in Eq. ^ are of order U, the final analytical result is 
correct to order U'^. We make an ansatz for p* similar 
toEq. (d]). 



Mc = 



Pc + ap^ 



1 



apl log 
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Eq. ^ shows that the "folding in" of F'pp^'^ for large 
frequencies iujm +ioJn2 gives a small contribution to 
This implies a reduced effective band width for the second 
order term, described by a different multiplying factor a 
in the logarithm. Identifying with the analytical result 
correct to order U'^, the factor a ss 0.10 is obtained. 
Eq. (fTO|) is then also correct to order U^. 
Fig. [1] shows results obtained by performing the calcu- 
lations in Eq. Q numerically using results for F'pp' up 
to first or second order in U and x'^ with S(iw„) = 0. 
The figure also shows the analytical result in Eq. (|T0|) . 
The second order result is clearly larger than the first 
order result. For pc ^ 0.5, Eq. pU|) describes the second 
order calculation rather well, while for larger pc correc- 
tions to the analytic result make p* still larger compared 
to Eq. (fT0| . The second order contribution is reduced 
by retardation effects, like the first order contribution. 
However, the reduction is substantially less efficient, as 
described by the small factor a. 

The Morel-Anderson theory [Eq. ([T])] makes two main 
predictions: (i) as D is increased for fixed pc and Wph, Pc 
goes to zero and (ii) as pc is increased for fixed D/wph, 
p* saturates at the value l/log(I?/a;ph). (i) remains true 
when the second order contribution F'pp^'^ is taken into 
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account, but the numerical result in Fig.[T]shows that (ii) 
is violated, i.e, fi* does not saturate as U is increased. 

We also conduct a complimentary analysis and extract 
results for //* from the spectral gap Agp of the supercon- 
ductor at r = 0. The approach is similar to the original 
work by Morel and Anderson Q , which included only the 
first order term in U . We work on the imaginary axis in 
the limit T — > 0. Starting point is the self-consistency 
equation for the off-diagonal self-energy, 



0.25 



S2l(iw„) 



(11) 



where the kernel K{iujn, iuJm) includes the attraction me- 
diated by the phonons and the repulsion to order C/^, 



A 
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Po 1 



+ U - U^U{iuJn+iuJm)- 



(12) 

and 



G2i(iuJm) is the offdiagonal Green's function [6|, 
a semi-elliptic DOS po{e) = VW^^/{2TTt'^) with W = 
At = 2D is used. The effect of the diagonal self-energy 
Ell is taken into account in the analytical calculations 
by a factor Z = 1 — E'j^j^(O), which is taken finite only for 
\oJn\ < ujph [see Eq. J^]. 

The self-consistency equation pT|) can be solved nu- 
merically by iteration. For an analytical solution, we 
need to make some approximations. At half filling, we 
use for the Green's function for |a;„| < Wph, 



G2l(lW„) 



1 



(13) 



for Wph < \ujn\ < D, G2i(iwn) ~ -E2i(ia;„)/(t|a;„|), and 
for \ujn\ > D, G2i{iuJn) — 0. A suitable ansatz for the 
off-diagonal self-energy is 



S2i(«w„) = A3 + A2/(iw„po) + 



Ai 



1 



(14) 



For /(.t) we use Eq. except that the numerical coef- 
ficients are modified for the semi-elliptic DOS. By com- 
paring the w-dependence with the full numerical solution 
we find reasonable accuracy for this functional form. We 
have to solve for the three parameters Ai, A2, and A3 
by evaluating the self-consistency equation at suitable 
values of ito. The general case is algebraically quite in- 
volved. Here, we only treat the first and purely second 
order cases explicitly to see the major effects. 

For the first order case, we set A2 = and omit the 
U^-term in Eq. IT^ . We use the conditions E2i(0) = 
Ai, and E2i(i-D) ~ A3, and assume A^ ^ Wph <C D. 
With the usual approximations we find a solution for the 
spectral gap Agp = Ai/Z of the form (23. 24 1. 



isp = ciWph exp 



ZC2 



A-m*(1 + C3A)-m:,i 



(15) 




FIG. 2: (Color online) The spectral gap Asp as calculated 
from the numerical solution of Eq. (Illf) with different kernels 
as a function of /ic for A = 0.5 and D/tUph = 80 in comparison 
with the corresponding analytical results based on Eq. (|15p . 



The result for /x* is given in Eq. ([T|) and fi* ^ =0. 

In the situation when only the C/^-term is included we 
set A3 = and omit the constant [/-term in Eq. ((T^). 
To determine the parameters Ai and A2, we use the 
following two conditions: E2i(0) = Ai, and 'Ej2i{'iD) + 
T,2i{—iD) ~ 2A2/(i/2). The calculation again yields a 
result of the form (jlSp . however, now with 



1 + log (tr) i + «/^?i°gte) 



9 4 



(16) 

As before < a2 < 1 accounts for the less effective re- 
tardation effects. In addition a higher order term /i* i 
appears, which was absent in the first order calculation. 
Such terms can account for the discrepancy between nu- 
merical and analytical results in Fig. [H where the nu- 
merical result does not saturate as function of /Xc. In our 
approximation the coefficient is 7 ~ 0.8 for D ^ Wph. 
In Fig. [51 we show the /ic-dependence of the numerical 
solution of Eq. pTjl for A = 0.5. It is compared with the 
analytical result in Eq. (fTS)) with the respective results 
for /!*, Eqs. (|l|10|16p . We use a as before and a2 ~ 0.2. 
ci — 1.7, C2 = 1.07 were determined by fitting to the 
numerical solution for = in the regime < A < 0.5. 
C3 = 0.8 was found to give a reasonable fit for the first 
order calculation in ^c- We omit the term //* ^ for the 
values of /Xc considered. The agreement between numer- 
ical and analytical results is quite good, which supports 
the earlier findings from the projection approach Eq. (j4|. 
One finds similar results as in Fig. [2]when calculating Tc. 
We now want to corroborate our findings with DMFT 
and analyze the impact of increasing /Xc on supercon- 
ductivity similar to Fig. [2l For this purpose we need 
to include additional effects. We extended our previous 
DMFT and ME perturbation theory (PT) Q to the case 
of finite U . In the PT for En and E21 we include direct 
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FIG. 3; (Color online) DMFT result for the spectral gap 
Asp ~ 2:E°*'(0) as a function of /ic for D — 2, constant A ~ 1 
according to the second order result for in Eq. (|17p and 
Eq. (|19|l . and lvq ~ 0.05 in comparison with the PT. We also 
included the result of the analytic formula in Eq. (|15p with 
renormalized parameters ojph — ojq, Z, calculated from PT, 
from Eq. (fTO)) and Ci as in Fig. [2l 

terms in U up to second order, for instance, the terms 
described in Eq. ([TT|). In addition, we have to take into 
account the fact that the electron-phonon vertex F^'^p) is 
renormalized by the Coulomb repulsion. We introduce 
the quantity r[^^\ which contains g and all corrections 
from the [/-term. For a weak frequency dependence up 
to the small phonon scale, we can define a renormalized 
coupling g'' = f[JP^(0,0). For g^ we use the result up to 
second order, 

— = 1 - aiHc + (a? - a,2)nl, (17) 
9 

where ai = -n(0)/po and a2 = ^ X)/, n(*'^fc)G'(icjfc)^. 
This was found to give a good description up to [/ ^ D 
[2^. Thus for the electron phonon part we use 

m 

(18) 

where the phonon propagator D{iuJm) is taken from the 
DMFT calculation [6|. The effective A is defined by 

00 

A = 2po[.9l'/d^^, (19) 


where the phonon spectral function p^{uj) includes self- 
energy corrections due to U. We define the renormalized 
phonon energy scale Wph = Wq by the peak position of 
p^{ijj). Through the condition A ~ 1 and Wph = ~ 
0.05 a set of bare model parameters (g,cJo,Mc) is deter- 
mined by the DMFT calculations, for which we can com- 
pare the PT with DMFT. The results are shown in Fig.H 



We find good agreement of the DMFT result with PT 
and the analytic formula, Eq. p^ . up to fic ~ 0.4 — 0.5. 
This demonstrates that (i) the electron-phonon vertex 
correction according to Eq. p7|) is suitable, (ii) that the 
higher order form for the Coulomb pseudopotential in Eq. 
PU| captures correctly the results of the PT and the full 
DMFT calculation, and (iii) that the effective parameter 
description is appropriate. Therefore, this validates the 
previous analysis in a more complete calculation and it 
corroborates our findings for ji* up to intermediate values 
of fic- For larger values of fic, we find that Agp in the PT 
calculations exceeds the DMFT result, where Agp 0. 
Then both Eq. dUl) for 5"" and Eq. ^ for p* start to 
underestimate the reduction effect. 

In conclusion, we have calculated corrections to the 
standard result for the Coulomb pseudopotential /i* . We 
provide numerical calculations of fi* and Asp as well as 
approximate analytical calculations. Whilst retardation 
effects still occur they are not as efficient for the higher 
order result [Eq. PU]) ]. Thus, for systems with sizeable 
effective Coulomb repulsion we should expect larger val- 
ues for /i* than what is traditionally quoted. This can 
help to explain the observations in a number of materials, 
such as Li, and is in line with findings for A from DFT 
calculations [l^ . 
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